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Abstract 

We show that the resummation of large radiative corrections in QCD 
processes can be performed both in covariant gauge and in axial gauge. We 
extend the resummation technique to inclusive processes, concentrating on 
deeply inelastic scattering, Drell-Yan production, and inclusive heavy meson 
decays in end-point regions. Leading and next-to-leading logarithms pro- 
duced by radiative corrections are organized into Sudakov factors, whose 
expressions are the same as those derived from the conventional Wilson-loop 
approach. The gauge invariance of the resummation in the end-point region 
is demonstrated explicitly. 
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I. INTRODUCTION 



One of the essential subjects in perturbative QCD (PQCD) is to organize 
large radiative corrections. For processes involving both a high and a low 
energy scale, large logarithms associated with ultraviolet, collinear and soft 
divergences are produced by loop integrals. Collinear and soft divergences 
may overlap to form double logarithms. These logarithms appear in products 
with coupling constant a s , and spoil the perturbative expansion. It has been 
found that radiative corrections in kinematic end-point regions of inclusive 
processes, such as deeply inelastic scattering, Drell-Yan production |^, [| 
and heavy meson decays B — > X u lv and B — > [|, |5j, [|, introduce double 
logarithms. Double logarithms also occur in exclusive processes, such as 
Landshoff scattering JF], hadron form factors ||, Compton scattering |J and 
heavy-to-light transitions B — > n(p) [0 and B — > D^*> at maximal recoil 



[ pnj] . In order to have a meaningful PQCD analysis of these processes, the 
important logarithms must be summed to all orders. 

The conventional approach to the summation of large radiative correc- 
tions in inclusive processes is based on the factorization of Wilson loops 
|T2] , 13, 13 1, which are formed by the classical trajectories of quarks involved 



in scattering. These Wilson loops collect important infrared logarithms, and 
are treated by renormalization group (RG) methods. An alternative ap- 
proach is the resummation technique developed in |TS) , which is more closely 
related to PQCD factorization theorems. This technique has been applied to 
exclusive processes mentioned above. Recently, it was extended to inclusive 
decays of heavy mesons || || , and results are found to be the same as those 
from the Wilson-loop formalism. 

In this paper we shall discuss the application of the resummation tech- 
nique to inclusive processes in a more rigorous and detailed manner, con- 
centrating on deeply inelastic scattering and Drell-Yan production. The 
equivalence of the resummation technique and the Wilson-loop formalism 
is justified by demonstrating that they lead to Sudakov factors which are 
equal up to next-to-leading logarithms. We shll show that a resummation 
can be performed both in covariant gauge d ■ A = and in axial gauge 
n ■ A = 0, and that the resummation is independent of the gauge vector n. 
This is a subject which cannot be addressed in the Wilson-loop approach, 
since the vector n has been fixed in the direction of classical trajectories of 
quarks. The gauge dependence of the resummation was not discussed either 
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in |TJ|], where a vector n with = was chosen. 

We show how to perform a resummation in covariant gauge and in axial 
gauge in Sec. II, taking a jet subprocess as an example. We then apply the 
technique to deeply inelastic scattering in Sec. III. The factorization of a 
relevant structure function into a distribution function, a jet function, a soft 
function and a hard scattering amplitude is formulated. Leading and next-to- 
leading logarithms in these convolution factors are organized into a Sudakov 
factor. It is justified explicitly that the dependence of the Sudakov factor 
on the gauge vector n cancels among the distribution, soft and jet functions. 
Drell-Yan production and inclusive heavy meson decays are studied in Sec. 
IV and Sec. V, respectively. We conclude that the resummation results for 
the above processes are the same as those from the Wilson-loop approach. 
Sec. VI is the dicussion. Some details of calculation are presented in the 
Appendix. 



II. COVARINT GAUGE AND AXIAL GAUGE 

The resummation technique was originally developed in axial gauge n-A = 
[U| . In this section we shall show that it can also be performed in covariant 
gauge d ■ A = 0, and leads to the same results. The physical meaning of the 
gauge vector n, and the concept of the resummation are in fact understood 
easilier in covariant gauge. 

We take as an example a jet subprocess defined by the matrix element 



J(p, n) = (0|Pexp 



DC 



ig / dzn ■ A(nz) 



n 



?(0)|P> , (1) 



where q is a light quark with momentum p. The abelian case of this sub- 
process has been discussed in flTj|. Our formulation presented below is more 
general and also suitable for the nonabelian case. The path-ordered expo- 
nential in Eq. (P is the consequence of the factorization of collinear gluons, 
possessing momenta parallel to p, from a full process. This factorization will 
become transparent in Sec. III. For convenience, we assume that p has a 
large light-cone component p + , and all its other components vanish. A gen- 
eral diagram of J is shown in Fig. 1(a), where the path-ordered exponential 



3 



is represented by an eikonal line along the vector n with collinear gluons 
attaching to it. The eikonal line behaves like a scalar particle carrying a 
color source, whose propagator is written as l/(n • /), I being the momentum 
flowing through the eikonal line. The Feynman rule for an eikonal vertex is 
given by n^T 0- , T a being a color matrix. All these rules can be derived from 
the definition of J in Eq. ([I]). To ensure that the collected gluons lead only 
to the collinear divergences from the region with I parallel to p, n does not 
lie on the light cone. 

It is easy to find that J contains double logarithms from the overlap of 
collinear and soft divergences by investigating the lowest-order diagrams in 
Figs. 1(b), the self-energy correction to the quark q, and 1(c), the vertex 
correction with a gluon connecting the quark line and the eikonal line. In 
covariant gauge both Figs. 1(b) and 1(c) produce double logarithms. In 
axial gauge n • A = 0, however, the path-ordered exponential reduces to an 
identity, and thus Fig. 1(c) does not exist. In this case n goes into the gluon 
propagator, (—i/l 2 )N tJ-u , with 



N^=g^ + n 2 —— 2, 2 



and only Fig. 1(b) leads to double logarithms. The vanishing of Fig. 1(c) in 
axial gauge can also be confirmed by the fact that the attachment of a gluon 
to an eikonal vertex gives a factor N^ v n v = 0. 

We show how to resum the double logarithms in J in covariant gauge. 
The essential step in the resummation procedures is to derive a differential 
equation p + dJ/dp + = CJ |5], ||, f[0f| , where the coefficient function C contains 
only single logarithms, and can be treated by RG methods. Note that the 
path-ordered exponential is scale-invariant in n, and thus J must depend on 
p and n through the ratio (p ■ n) 2 /n 2 . Hence, the differential operator d/dp + 
can be replaced by d/dn using a chain rule, 

, d n 2 d , . 

i r——j= Ua _j (3) 

dp + v ■ n dn a 

v = (1, 0, Ot) being a vector defined by p = p + v. 

Equation @ simplifies the analysis tremendously, because n appears only 
in the Feynman rules for the eikonal line, while p may flow through the whole 
diagram in Fig. 1(a). The differentiation of each eikonal vertex and of the 
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associated eikonal propagator with respect to n, 



v ■ n dn n n ■ I 



n 2 d n 




) 



n ■ I 



1 



(4) 



leads to the derivative p + dJ/dp + described by Fig. 2. The summation in 
Fig. 2 is performed over different attachments of the symbol +, which denotes 
the new vertex n M defined by the second expression in Eq. (|J). 

If the loop momentum / is parallel to p, the factor v ■ I vanishes, and n M is 
proportional to v^. When this n M is contracted with a vertex in J, in which 
all momenta are mainly parallel to p, it can be shown that the contribution to 
p + dJ/dp + from this collinear region diminishes. Hence, the leading regions 
of I are soft and hard. 

According to this observation, we investigate some two-loop examples ex- 
hibited in Fig. 3(a). If the loop momentum flowing through the new vertex 
is soft but the other one is not, only the first diagram is important, giving 
large single logarithms. In this soft region the subdiagram containing the 
new vertex can be factorized using the eikonal approximation as shown in 
Fig. 3(b), where the symbol ® represents a convoluting relation. The sub- 
diagram is absorbed into a function K, and the remainder assigned to J, 
both being 0(a s ) contributions. The function K is part of the coefficient C 
introduced before. If both the loop momenta are soft, the four diagrams in 
Fig. 3(a) are equally important. The subdiagrams are factorized according 
to Fig. 3(c), which give 0{a 2 s ) contributions to K. The remainder is then 
the leading diagram of J. 

If the loop momentum flowing through the new vertex is hard and the 
other one is not, only the second diagram in Fig. 3(a) is important. In this 
region the subdiagram containing the new vertex is factorized as shown in 
Fig. 3(d). The right-hand side of the dashed line is absorbed into a function 
G as an 0(a s ) contribution, which is another part of the coefficient C . The 
left-hand side is identified to be the 0(a s ) diagram of J. If both the loop 
momenta are hard, all the diagrams in Fig. 3(a) are absorbed into G, giving 
0(a 2 ) contributions. 

Extending the above analysis to all orders, we derive the differential equa- 
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tion 

P+ dp+ J= - ft '( m // i ' as (/ i )) + G {P +U I ^ a s{^)) J, (5) 

with the functions K and G described by the corresponding diagrams in 
Fig. 4(a). In Eq. @ is a factorization scale, and the gauge factor v = 
\J [y ■ n) 2 /\n 2 \ in G comes from the reexpression J {p ■ n) 2 /\n 2 \ = p + u. A 
gluon mass m is introduced to regulate the infrared divergences in K. We 
have made explicit the fact that K and G depend on a single infrared scale 
m and a single ultraviolet scale p + , respectively. 

We consider the 0(a s ) contribution to K from Fig. 4(b), which is written 

as 

d 4 ~ e l n„ 
(27r) 4 - e n-//2 

f d^~ e l n 2 
= -^^1 (27r) 4- e (n . /)2(/2 _ m2) " SK > ( 6 ) 

where Cp = 4/3 is a color factor and SK an additive counterterm. The 0(a s ) 
contribution to G from Fig. 4(c), where the soft subtraction ensures a hard 
loop momentum flow, is given by 

G = -ig C Ff i / , , 4 _ f „ z -n 7- , A2 7^ - — -7 - <*<3 , (7) 



• 2/, e [ d?*l g* 1 " v„ 



(27r) 4 " e n-/ Z 2 \(P + 2 vi, 
being an additive counterterm. By evaluating Eqs. (^) and (0) di- 
rectly, one finds that K and G contain only single logarithms ln(m//x) and 
ln(p + z//yu), respectively, as claimed before. Organizing these single loga- 
rithms using RG methods, and then solving Eq. @, we resum the double 
logarithms in J. These procedures will be implemented in a detailed way in 
Sec. III. 

We then review how to perform a resummation in axial gauge n ■ A — 0, 
and compare the outcomes with those obtained in covariant gauge. Since 
gluons are detatched from the eikonal line in axial gauge, we concentrate on 
diagrams of the type Fig. 1(b). In this case the scale invariance of J in n is 
manifest in the gluon propagator N^ v . The chain rule in Eq. @ then still 
holds, but now d/dn operates on N^, giving 

2 A 2 



v ■ n ° dn a v ■ nn ■ I 

= v a (NT + N au F) . (8) 
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The momentum V 1 (l u ) appearing at one end of the differentiated gluon line is 
contracted with a vertex the gluon attaches, and the vertex is replaced by the 
new one v a , defined by the second expression in Eq. (|8|). The contraction of 
I 11 hints the application of the Ward identity Summing all the diagrams with 
different differentiated gluon propagators, and employing the Ward identity, 
one finds that the new vertex moves to the right most end of the g-quark 
line. We then obtain the derivative p + dJ/dp + described by Fig. 5(a), where 
the square represents the new vertex v a . 

We argue that the collinear region of the loop momentum I is not im- 
portant because of the factor l/(n • I) in v a with nonvanishing n 2 . Hence, 
the leading regions of I are soft and hard, in which the subdiagram con- 
taining the square can be factorized into a function K and a function G, 
respectively, following the same reasoning as for Fig. 3. We then derive a 
differential equation similar to (^) with the general diagrams of K and G 
displayed in Fig. 5(b). The 0(a s ) contributions to 1C from Fig. 5(c) and to 
Q from Fig. 5(d) are given by 



K = ig 2 Cp[Ji t J 



5K 



(27r) 4 - e (vl)(l 2 -m 2 ) 
f d 4 ~ e l n 2 

J (2n)*~< (n • W ~ m 2 ) " * K ' (9) 

d*~n „ ( J±J_ 

(2tt) 4 -^ Up + If 



G = ig 2 C F ^j 7 — —v fl [ 7 ^-J^j v --^.)—-SG. (10) 




Similarly, a gluon mass m is introduced to serve as an infrared cutoff of the 
integral for K. The second expression in Eq. (H) is derived in the Appendix: 
The first term g^ u in does not contribute. The contribution from the 
third term —n u l^/(n ■ I) cancels that from the fourth term •n 2 \H v /(n ■ I) 2 . 
Hence, only the second term —nH v j {n ■ I) remains. 

It is observed that Eq. (|) in covariant gauge is equal to @ in axial 
gauge. It is also straightforward to justify that Eq. (|7|) is equal to flrO|). 
These equalities generalize to higher loops, and thus the differential equation 
(H) in covariant gauge is exactly the same as that in axial gauge. We conclude 
that the resummation can be performed in both covariant and axial gauges, 
and the same results are obtained. In the following analyses of inclusive 
processes we shall work only in axial gauge. 
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III. DEEPLY INELASTIC SCATTERING 



In this section we discuss the resummation of large radiative corrections 
in deeply inelastic scattering, 

£(k) + h{p) — > £(k') + X , (11) 

where k (k 1 ) is the initial (final) momentum of the lepton £, and p is the 
momentum of the initial hadron h. We choose a center-of-mass frame of the 
incoming hadron and lepton, such that p and k have the light-cone compo- 
nents p = ip + , 0, Ot) and k = (0, k~, Or) with p + = k~ . The Bjorken variable 
is defined by 

-f^^, (12) 

q = k — k' being the momentum tranfer. We shall derive the factorization 
formula for deeply inelastic scattering, which is appropriate in the end-point 
region x — > 1. Radiative corrections to the convolution factors involved in 
the factorization formula that produce double logarithms will be organized 
using the resummation technique in axial gauge. 

III. 1 Factorization 

Fig. 6(a) describes the scattering of a valence quark in the hadron h with 
the fractional momentum £p + Pt by an off-shell photon with the momentum 
q. Pt is a small transverse momentum, which will serve as the infrared 
cutoff of loop integrals below. The outgoing quark carries the momentum 
p' = £p + Pt + q — {£,P + ~ k' + , k~ — k'~, px — k T ), where the components 
satisfy the on-shell condition k' 2 = 0. Including the requirement p' 2 > 0, we 
have the constraint x < £ < 1. 

We consider radiative corrections to Fig. 6(a) in axial gauge. Some 0(a s ) 
diagrams are displayed in Figs. 6(b)-(e). Fig. 6(b), the self-energy correction 
to the valence quark, and Fig. 6(c), the loop correction with a real gluon 
connecting the two valence quarks, contain both collinear divergences from 
the loop momentum / parallel to p and soft divergences from small I. Since 
the soft divergences cancel asymptotically as shown later, the double loga- 
rithmic corrections are mainly collinear. Therefore, Figs. 6(b) and 6(c) can 
be absorbed into a distribution function associated with the initial hadron. 

In the end-point region with k l+ — > p + and k'~ — > (ie., x — > £ — > 1), the 
outgoing quark, possessing a large minus component p'~ ~ k~ but a small 
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invariant mass p' 2 ~ (£ — x)Q 2 , produces a jet of collinear particles. The 
self-energy ccorrection to the outgoing quark in Fig. 6(d) then gives both 
collinear divergences from I parallel to p' and soft divergences. Because of 
the asymptotic cancellation of the soft divergences, these double logarithms 
are mainly collinear, and Fig. 6(d) can be absorbed into a jet function J. 
Note that this jet function does not exist in the conventional factorization 
formulas for deeply inelastic scattering, which are appopriate for x not close 
to 1. This is the reason the logarithmic corrections are stronger in the end- 
point region. 

The absorption of the vertex-correction diagram in Fig. 6(e) is more del- 
icate. Fig. 6(e) contains collinear divergences from / parallel to p and soft 
divergences. Moreover, it contains additional collinear divergences from I 
parallel to p' as x — > 1. We treat this diagram in different ways in these 
three leading regions. For / parallel to p, we replace the outgoing quark by 
an eikonal line in the direction n, and factorize the gluon into from the 
full scattering amplitude. It is then clear that such an eikonal line, collect- 
ing collinear gluons, is necessary for the factorization of vertex-correction 
diagrams, and has been associated with the jet subprocess in Eq. ([J). This 
arbitrary vector n with nonvanishing n 2 does not introduce extra collinear 
divergences, but soft divergences through the eikonal propagator l/(n • I). 
These soft divergences, differing from the original ones which should be 
grouped by the classical trajectory of the outgoing quark in the direction 
p', can be removed accorrding to Fig. 7(a). The soft subtraction is obtained 
by further replacing the valence quark by an eikonal line in the direction 

u=p/p+ = (l,0,Or). 

Similarly, for I parallel to p', we replace the valence quark by an eikonal 
line also along n, and absorb the gluon into J. The additional soft diver- 
gences are subtracted following Fig. 7(b), where v' = (0, 1, Ot) is a vector 
along p' in the x — > 1 region. If I is soft, all quarks are replaced by their 
classical trajectories, ie., by an eikonal line along v for the valence quark and 
by an eikonal line along v' for the outgoing quark, as shown in Fig. 7(c). 
This diagram, together with the collinear subtractions which remove the di- 
vergences from / parallel to v and to v ', are factorized into a soft function U. 
Fortunately, in axial gauge n • A = all the diagrams containing eikonal lines 
in the direction n do not contribute, and only the first diagram in Fig. 7(c) 
is left. At last, when the loop momenta in Figs. 6(b)-(e) are hard, the gluons 
are grouped into a hard scattering amplitude H. 
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Based on the above reasoning, we write down the factorization formula 
for a structure function F(x, Q 2 ) associated with deeply inelastic scattering, 



F(x,Q 2 ) = J\ J dM{^h^)J{p',b^)U{b^)H{^Q^) , (13) 

graphically described by Fig. 8, where b is the conjugate variable of pt- The 
reason Eq. (|i~3"D is formulated in b, instead of pt, space has been given in 
[^, TUfl . Briefly speaking, a formulation in b space facilitates the factorization 
of higher-order corrections. The variable b is introduced by a Fourier function 
e llj " h associated with the integrand for a higher-order correction whenever 
the loop momentum I flows through <ft or J. We have assumed that all the 
convolution factors in Eq. ( p~3"D depend on a single b for simplicity. The 
relevant scales each convolution factor contains have been shown explicitly. 
<f) and J, involving both a large scale p + (or p'~) and a small scale 1/b, 
possess double logarithms. U and H, depending only on one scale 1/b and 
Q, respectively, contain single logarithms. 

III. 2 The distribution function </> 

We demonstrate how to resum the double logarithms in the distribution 
function (p. For a similar reason, <fr depends on the ratio (p ■ n) 2 /n 2 , and 
thus the chain rule in Eq. (|3|) holds. Applying d/dn to a gluon propagator 
and summing all the diagrams with different differentiated gluon lines, we 
obtain the derivative p + d(j)/dp + described by Fig. 9(a), where the coefficient 
2 counts the outer most ends of the two valence quark lines, and the square 
represents the new vertex v a defined by Eq. (H). 

Similarly, in the leading soft and ultraviolet regions the subdiagram con- 
taining the square is factorized into a function K and a function G, respec- 
tively. The differential equation is then expressed as 



P 



d 
dp^ 



K(bfj.,a 8 (n)) + G(p + u/n,a s (n)) 0, (14) 



with the general diagrams of K and G shown in Fig. 9(b). We keep the gauge 
factor v = J {v ■ n) 2 /\n 2 \ in G in order to trace the gauge dependence of the 
resummation. It has been shown in Eq. flliD that K depends on a single 
small scale 1/b and G depends on a single large scale p + . 
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The 0(a s ) contributions to K from Fig. 9(c) and to G from Fig. 9(d) are 
given by 



K 

G 



ig 2 C F /j e 
ig 2 C F /j l t 



(2vr) 4 - 

d 4 -n 



- + 2inS(l 2 )e llT - h 



V J^N^-5K, (15) 
v ■ I 



/ 



(2vr) 4 -^ \(p + l) 2lu v-l) I 2 



" 1 - SG 



(16) 



SK and 8G being the corresponding additive counterterms. In Eq. (|T5| ) the 
second term, corresponding to the second diagram in Fig. 9(c), contains the 
factor e dT ' b , such that K is free of infrared poles. That is, 1/6 serves as an 
infrared cutoff of the loop integral as stated before. Note that K and thus 
the soft divergences vanish in the asymptotic region 6^0. The second term 
in Eq. ( [To]) is the soft subtraction associated with G. 

A straightforward calculation of Eqs. (|15|) and (|l^) in the modified min- 
imal subtraction (MS) scheme leads to 



K(bfj,,a s ) 
G(p + u/fi,a s ) 



— —Cp ln(6//) . 

7T 

-^C F ln^ 

7T jl 



(17) 
(18) 



and 8K = —5G. Details of the calculation refer to the Appendix. To perform 
the integral for G, we have chosen a space-like gauge vector n with n 2 < 0, 
which avoids imaginary contributions. In Eqs. (0) and (pl|) we keep only 
the large logarithmic terms, and neglect other constants of order unity. This 
is reasonable, because the factorization of K and G in Eq. ([14]) holds only 
in the leading regions characterized by large logarithms. It is found that G 
depends on the large scale p + through the ratio (p ■ n) 2 /n 2 as expected. 

Since K and G contain only single soft and ultraviolet logarithms, respec- 
tively, they can be treated by RG methods: 



li-^-K = -X K 
an 



d 

-li— G 

djj, 



(19) 



The anomalous dimension of K, Xk 

by EE 



fid5K/dfj,, is given, up to two loops, 



with n/ the number of quark flavors, and Ca = 3 a color factor. In solving 
Eq. ([H^), we allow the scale /i to evolve to the infrared cutoff 1/6 in K and 
to p + in G. The RG solution of K + G is then written as 

if(6/i, a s (//)) + G(p + u/fi, a s (n)) 



ib p 

). . fP' dp,. . ... 
-C F \nu- / — Ajr(a a (/x)) . 

Jl/6 « 



Substituting Eq. 

<f>(£,P,b } fJ,) = exp 



a s (l/6)) + G(z/, a s (p + )) - / ^A^(« s (/i)) 

Ji/fe /i 

/6 /U 

1|) into (|14j) , we derive 
Q dp 



(21) 



1/6 p -1 



— X K {a s {p)) H C F lnz/ 

1/6 /i 7T / 



x0(f,M, (22) 

We have set the lower limit of the variable p + to 1/6, and the upper limit 
to Q, which is of the same order as p+ ax = Q/^/2. The initial condition 
</>(£, 6, /x) contains single ultraviolet logarithms from self-energy corrections to 
the valence quarks, which can be further organized through the RG equation, 



with 



£>0(£,M = -2A g 0(£,M 



d d 



(23) 



(24) 



and \ q = —OL s j 7r the quark anomalous dimension in axial gauge. The solution 
to Eq. is given by 



</>{^,b,fj) = exp 



-2 r ^\ q (a s m 

Jib U 



0(£,6,l/6). 



(25) 



^ dp 
'1/6 /2 

Combining Eqs. (p^) and (p5[), and employing the variable changes p + = 
in the first integral and p 2 = yQ 2 in the second integral, we resum the 

large radiative corrections in the distribution function cf) up to next-to-leading 

logarithms: 



<!>(£, P,b,p) 



exp 



1 dy 

i/(6Q) y 



— Ajf («,(//)) + 2 C F In i/ 

1/6 2 jjL z 7T , 



l/(60)a y 



0(^,6,1/6). 



(26) 
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III. 3 The jet function J 

We now discuss the resummation in the case of the jet function J. There 
exists a similar chain rule for J, 

/ d n 2 . d T . . 

p V = -^ J ' (27) 

Employing the same procedures as in Sec. II, we obtain the derivative 
p'~dJ/dp'~ described by Fig. 10(a), where the coefficient 2, again, comes 
from the two outer most ends of the outgoing quark line, and the square 
represents the new vertex 

v' a = -7^7 • (28) 
v' ■ nn-l 

The subdiagram containing the new vertex is factorized in the leading soft 
and ultraviolet regions into a function K' and a function G' , respectively, 
whose general diagrams of K' and G' are shown in Fig. 10(b). 

The resummation of large radiative corrections in J is subtler. We argue 
that the resummation for results in Sudakov suppression in the large b 
region, and thus b can be set to zero approximately in the analysis of J. 
Another reason for this approximation is that the function K' can be made 
free of infrared poles without introducing a b. This is obvious from the 0(a s ) 
expression of K' from Fig. 10(c), 



5(p' 2 )K' = ig 2 C F ^ J 



(27r) 4 " e v' ■ I 



S(p' 2 ) 
I 2 



+ 2m5(l 2 )5((p' -iy 



S(p' 2 )5K' , (29) 

where the 5 functions correspond to the final-state cut, and SK' is an additive 
counterterm. It is easy to observe that p'~ in the second term regulates the 
infrared divergence of the first integral. In order to extract the on-shell pole 
of the outgoing quark, we apply the operator 

/+ 

5{p' 2 ) I dp' 2 = 5{p' 2 ){2p'-) [ Puia *dp' + (30) 



max 



to the second term of the loop integral. In the case of K' we choose p 
p'~ . After a tedious calculation as exhibited in the Appendix, we have 

K' = -^C F \n^ (31) 

7T p~ 
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with v' — J (v' ■ n) 2 /\n 2 \ the gauge factor. 

The 0(a s ) expression of G' from Fig. 10(d) is given by 



G' 



-ig 2 C F fx e 



Cf In 

7T 



d A ~n „, 

( 27r )4-e U M 
1 p+ N 



7* 



/ 



z/ ^Jp + p'~ 



(p'-l) 2 v'-l 
SG 1 , 



(32) 



where is an additive counterterm. In the above formula the loop integral 
comes from the first two diagrams in Fig. 10(d), and the second term from the 
third and fourth diagrams. We obtain the second term by simply choosing 
Pmax = P + ~ Q m Eq. (|30"D, which is greater than p^" ax = p'~ for the soft 
function K', and by choosing yfq 1 ~ y/p + p'~ , which is the typical scale of the 
third diagram, as the infrared cutoff of the integral for the fourth diagram. 
A similar calculation leads to 



G' 



7T H 



Hence, the general differential equation for J is written as 
d 



P 



dp'' 



-J 



K'{v'n/p' , a s (fi)) + G'(^p + p'-//j,, a 8 (ji)) 



(33) 



J. (34) 



The functions K' and G' for J possess the same functional forms as K and 
G for (J), respectively, and the counterterms have the relations SK = SK' and 
5G = SG'. We organize the single logarithms in K' and G' to give 



K'(v'n/p'-, a e (ji)) + G(^/p+p'-/ii, a s (/x)) 



ot s {p' 



P~) n , / /•v 7 ^ dp, 
C F \nu - / — X K (a s {fx)) . 

71 Jp'- fl 



(35) 



Solving Eq. 



we obtain 



J(p', b, //) = exp 



Q dp' 



oi s {p'~ 



7T 



1/6 p' 



-Cp In v' 



— \ K {a s {\i)) 



J(b,fi) . 



(36) 
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The initial condition J(b, //) obeys a RG equation similar to Eq. (|23|), and 
is thus given by 



J(b, //) = exp 



f * ^L-\ q (a s (jl))\ Jib, y/Q/b) , (37) 

JQ/b /i 



which sums the single ultraviolet logarithms arising from self-energy correc- 
tions to the outgoing quark. We have chosen the invariant mass of the jet 
p' 2 ~ p'~p' + ~ Q/b as the lower bound of jl 2 in the above integral. Com- 
bining Eqs. (p6|) and ([37]), we resum the large radiative corrections in the jet 
function J up to next-to-leading logarithms: 

J(p', b, /i) = exp - / ^ ( f[ Q ] ^AkKG/)) + 2^1c F In z/) 




where the variable changes similar to those for <fi have been employed. Note 
that we have taken into account four scales in the resummation procedures 
for J, which have the ordering Q > ^/p + p'~ > p'~~ > 1/b > Aqcd- It is 
easy to observe that the integral generating double logarithms vanishes in 
the y — > 1, ie., p' ^> 1/6, region, where the jet is energetic. 

III. 4 The soft function U and the hard scattering H 

We have organized the single ultraviolet logarithms in the distribution 
function and in the jet function J into the exponentials in terms of the 
anomalous dimension X q . To sum up next-to- leading logarithms, we must 
include the single soft logarithms in U. We shall show that the gauge factors 
v and v' in and J, respectively, are canceled by those in U. 

The RG equation for U is written as 

VU(b,n) = -\vV(b,n) , (39) 

where the anomalous dimension is derived from the first diagram of 
Fig. 7(c). The one-loop bare soft function is given by 



U® = -2ig 2 C F ^f 



d 4 ~ e l u/JV" 
(2tt) 4 - £ vlv'- ll 2 



a 1 

-2— lmW)-. (40) 

7r e 
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The coefficient 2 in the first expression corresponds to the two vertex-correction 
diagrams, and those constants of order unity irrelevant to the gauge depen- 
dence have been neglected in the second expression. The gauge dependent 
factor lnz/ (lnz/) in Eq. (|40|) comes from the term nH v /n ■ I {n u l^/n ■ I) in 
N^ v . 

It is then trivial to compute \u from the multiplicative renormalization 
constant Zy, defined by = ZjjXJ which is written as 



Ac/ 

We obtain the solution to Eq. 

U(b, n) = exp 



-2— C F ln(iV) . 

7T 



p-IQ dy 

—\u(a s (yQ)) 

i/(«j) y 



U(b, 1/6) . 



(41) 



(42) 



Since the structure function F is independent of fi, we have the RG 
equation for the hard scattering H 



VH((,Q,n) = (4\ 9 + \ u )H(t,Q,») 
which is solved to give 



(43) 



H(£,Q,/j) = exp 



dy 



v/q 2 y 

xH(£,Q,Q). 



dy 



im/q y 



^u{a s {yQ)) 



(44) 



In the above expressions for U and H necessary variable changes have been 
made. 

Substituting Eqs. Q3§), fl§3) and (g§ into QT3|), we derive 

F{x, Q 2 ) = £ di J d 2 b4>(Z, 6, 1/6) J(6, \[oJb)U{b, l/b)H{£, Q, Q) 

xexp[-S ms (Q,b)] , (45) 

where the Sudakov exponent Sbis is given by 

dy 



S 



DIS 



i/(6Q) y 
1 



yQ 2 da 2 

—X K (a s (fj,)) + \ q (at a (y/yQ)) 

l/b 2 /T 



+ 



dy 



i/(6Q) 2 y 



(46) 
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It is obvious that the gauge factors v from G and v' from K' have been 
canceled by those from U, and that the resummation technique developed 
above is gauge independent. This is an issue which cannot be addressed in the 
conventional Wilson-loop approach, since the vector n has been fixed in the 
direction v or v' in order to form a loop along the classical trajectories of the 
quarks. It is easy to observe that Eq. (f46|) is the same as the corresponding 
formulas derived from the Wilson-loop formalism with the identification of 
the anomalous dimension Xk = A in and Xk = r cusp in || . Note that the 
second integral involving X q in Eq. ( [46|) does not exist in [§, |], because the 
self-energy corrections to the valence quarks can not be taken into account 
by a Wilson (eikonal) line. The self-energy corrections to the outgoing quark, 
denoted by the second term in the first integral, also appear in || which 
were introduced, however, through an Altarelli-Parisi evolution equation. 



IV. Drell-Yan Production 

In this section we discuss the resummation for Drell-Yan production, 

h( Pl ) + h(p 2 )^£l(q)+X , (47) 

where p\ = (pf, 0, Ot) and p 2 = (0,^2 , 0^) with p± = p^ are the momenta of 
the incoming hadrons h in their center-of-mass frame, and q is the momentum 
of the lepton pair produced in the process. 

We start with the annihilation at the quark level as shown in Fig. 11(a), 
and then consider 0(a s ) corrections from Figs. ll(b)-(e). Figs. 11(b) and 
11(c) have been analyzed in Sec. Ill, and the conclusion is the same: They 
produce double logarithms with collinear ones dominant, and are thus ab- 
sorbed into distribution functions <p\ or 2 . The treatment of Figs. 11(d) and 
11(e) is similar to that of Fig. 6(e). If the loop momentum I is parallel to 
Pi, we replace the second quark by an eikonal line in the direction n, and 
factorize the gluon into <p\. Certainly, a soft subtraction with the first quark 
further replaced by an eikonal line in the direction v% = Pi/pf, similar to 
Fig. 7(a), is necessary. If I is parallel to p 2 , the diagram is factorized into 
02 following the above prescription. If / is soft, the first and second quarks 
are replaced by eikonal lines along v\ and v 2 = P2/P2, respectively, and 
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the diagram is factorized into a soft function U. A corresponding collinear 
subtraction similar to Fig. 7(c), should be included. 

In axial gauge n ■ A = contributions from all the above diagrams vanish, 
except for the one with both quarks replaced by eikonal lines along their 
classical trajectories. At last, when / is hard, Figs. ll(b)-(e) are absorbed into 
a hard scattering amplitude H. Hence, we arrive at the general factorization 
picture for the cross section of a Drell-Yan process as shown in Fig. 12, and 
the corresponding formula, 

,da(r,Q^ 



Q 2 C ' jq 2 — = J d£id£ 2 J ti 2 b0i(^i,pi,fo ) ^)0 2 (6,P2,&,A*) 

xU(b,fj,)H(Z u &,Q,ti), (48) 
where the variable r is defined by g 

O 2 

T = *-, (49) 
s 

with s = (pi + p 2 ) 2 and Q 2 = q 2 the invariant mass of the lepton pair. We 
have assumed that the first (second) quark carries a fractional momentum 
CiPl + Pit (£2^2 + P2t) with a small transverse momentum pxr (P2t), which 
regulates the infrared divergences of higher-order corrections. The cross sec- 
tion a depends only on a single b = b\ = b 2 , where b\ and b 2 are the con- 
jugate variables of pyr and p 2 T, respectively. The b dependence comes from 
the Fourier factor e llj " h , which is associated with a loop integrand whenever 
the loop momentum flows through <fti or <p 2 . 

We discuss the resummation of large radiative corrections in Drell-Yan 
production. We emphasize that double logarithms exist for arbitrary inter- 
mediate and large values of r, and no new double logarithms occur in the 
end-point region with r — > 1. This differs from the case of deeply inelastic 
scattering, where extra double logarithms appear at the end point. In order 
to compare the results for deeply inelastic scattering and for Drell-Yan pro- 
duction, we concentrate on the r — > 1 region, where both £1 and £2 approach 
unity. The resummation for the distribution functions 0i and <p 2 is the same 
as that in deeply inelastic scattering. We quote Eq. ( p6[) directly here: 

Jl/(bQ) V Wl/b 2 U Z 7T , 



>i{&,Pi,b,ii) = exp 



i/(6Q) y \Ji/b2 n 

v 2 /Q 2 dy 

—X q (a s (^/yQ)) 
i/(6Q)2 y 



&(£i,M/&) , (50) 
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for i — 1, 2, where the gauge factors z/j are defined by vi = y (u» • nj^/ |n* 

We then consider the soft function U . The lowest-order diagrams con- 
tributing to U are those with a gluon attaching two eikonal lines on the same 
side of the final-state cut. The diagrams with a real gluon attaching two 
eikonal lines on the different sides of the final-state cut are ultravioletly fi- 
nite due to the factor e llr ' b associated with the corresponding loop integrand. 
Hence, the anomalous dimension 



Af/ = —2— C F mfi^z^) 

7T 



(51) 



is similar to that for deeply inelastic scattering in Eq. ( |41| ) . The evolution of 
the soft function is then described by Eq. (|42|). Accordingly, the summation 
of the single ultraviolet logarithms in H leads to an evolution similar to 
Eq. ©. 

Combining the above results for each convolution factor, Eq. ( pE8D is re- 
expressed as 



Q 



2 da(r,Q 2 ) 



dQ* 



d£id& J d 2 bMZi,b,l/b)M&,b,l/b) 
xU(b, l/6)#(£i, Q, Q) exp[-S DY (g, b)} 



where the Sudakov exponent Soy is given by 



Si 



DY 



2 ^ 



dy fy*Q*dp,\ , / \ \ o 
" ' ^-A^(a s Gu)) + 2 



i/(6Q) y Ji/b* n 



l/(60)a y 



—\{u s {^/yQ)) 



(52) 



(53) 



Again, the gauge factors ^ from G have been canceled by those from U, such 
that the resummation technique is gauge independent. 



We consider the ratio A 

"i dy ' 
n/(bQ) y 



A 



exp 



2 ^ 



e- 5DY /(2e 

) 2 At.I 



2/Q ci/i 2 3 C F 

—X K (a a w) - o— a s{s/yQ) 



'v 2 Q 2 ^ 2 7r 

which is similar to the coefficient function Ajv in 



(54) 



A 



N 



exp 



dx — 



x 



N-l 



X 



—X K {a s {ii)) 

(1-X)2Q2 /i 2 



■~—a s (Vl - xQ) 

I 71 



(55) 
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with bQ corresponding to the large number N. is associated with the evo- 
lution of the iV-th moment of the differential cross section. Eq. ([54]) is exactly 
the resummation result of the jet function J in deeply inelastic scattering, 
excluding the gauge dependent factor. We conclude that the resummation 
technique is equivalent to the Wilson-loop approach in the study of deeply 
inelastic scattering and Drell-Yan production. 



V. INCLUSIVE HEAVY MESON DECAYS 

The technique developed above can be applied to inclusive heavy meson 
decays, such as B — > X s 7 and B — > X u lv, straightforwardly. The resum- 
mation of large radiative corrections in these processes at the end points of 
spectra |L7|, [18], |19| has been performed in || |[], where the specific gauge 
vectors n oc (1,1, Ot) and n oc (1,-1, Ot) were employed, respectively. In 
this section we shall discuss the case with a general n, and show that the 
results are gauge independent. 

The Wilson-loop formalism has been applied to the study of end-point 
singularities in inclusive B meson decays [|J based on the relation between 
Wilson lines and the heavy quark effective theory (HQET) G(J explored in 



2341 - The Wilson loop is absorbed into a heavy quark distribution function, 
whose evolution is governed by RG equations. We shall demonstrate that 
the resummation technique is equivalent to the Wilson-loop formalism in the 
study of inclusive heavy meson decays. 

We concentrate on the inclusive radiative decay B — > X s 7, which occurs 
through the transition b — > 57 described by an effective Hamiltonian |2~2[| . 
The basic factorization of this process at the quark level is similar to Fig. 6(a). 
The momentum of the b quark is P — k, P = [Mj v2)(l, 1, Ot) being the B 
meson momentum, and k being the momentum carried by light partons in 
the B meson, k has a plus component k + and small transverse components 
kj-. These momenta satisfy the on-shell conditions (P — k) 2 « M 2 and 
k 2 « 0. Here M and Mj, are the B meson mass and the b quark mass, 
respectively. The b quark decays into a real photon of momentum q and 
a s quark of momentum p', which is regarded as being light. Assuming 
that the nonvanishing component of q is q + , we have p' = P — k — q = 



20 



(P + -k+-q+,p-,-k T ). 

The absorption of infrared divergences is analyzed below. The self-energy 
correction to the b quark and the loop correction with a real gluon connecting 
the two b quarks as in Fig. 6(b) and 6(c), respectively, which contain only 
single soft logarithms, are grouped into a distribution function 0. In the 
end-point region with q + — > P + , the s quark has a large minus component 
p'~ = P but a small invariant mass p' 2 . The case is then similar to deeply 
inelastic scattering at the end point, and thus extra double logarithms appear. 
Therefore, the self-energy correction to the s quark in Fig. 6(d) is absorbed 
into a jet function J. 

The photon vertex correction gives both collinear logarithms from the 
loop momentum I parallel to p' in the end-point region and soft logarithms 
from vanishing /. If Z is parallel to p', we replace the b quark by an eikonal 
line in the direction n, and factorize the diagram into J. In axial gauge 
n ■ A = this diagram along with its soft subtraction do not contribute. If 
I is soft, the b and s quarks are replaced by eikonal lines along v = \/2P/M 
and v' = p' /p'~ , respectively ||. The replacement of a b quark by an eikonal 
line is consistent with the flavor symmetry in the HQET. We then absorb the 
gluon into a soft function U. Similarly, the collinear subtraction associated 
with this diagram does not contribute in axial gauge. 

Remaining important corrections from hard gluons are absorbed into a 
hard scattering amplitude H . At last, the factorization formula for the spec- 
trum of the decay B — > X s ^ in b space is written as 

dT r r 

^L = jdZj d 2 b<p(£, b, fi) J(p', b, fi)U(b, //)#(£, M, //) , (56) 

which is also described by Fig. 8. Here b is the conjugate variable of fcy, 
E 1 = q + 1 v^2 is the photon energy, £ = 1 — k + / P + is the momentum fraction, 
and /i is the factorization and renormalization scale. The resummation for 
J has been performed in Sec. II, and the result is the same as Eq. ( |38"D but 
with M substituted for Q. 

The single soft logarithms in <ft can be collected by eikonal lines in the 
direction v H. The self-energy correction to the b quark vanishes under this 
eikonal approximation. Hence, we consider only the diagram with a real 
gluon attaching the two eikonal lines, which leads to the loop integral for the 
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bare distribution function (fft\ 

^*W|^w 2rf( ' VT " b - (57) 

Evaluating Eq. (|57|), we obtain the pole term, and thus the anomalous di- 
mension of (f). 

Similarly, the bare soft function is written as 

u*» = I , d ':'! v »: lN z . (58) 

where the two vertex-correction diagrams have been included. The term 
—n v l^/ (n ■ I) in N^ u gives the gauge dependent factor In z/, which cancels 
that of J as in the case of deeply inelastic scattering. The contribution from 
the term nH v / (n ■ I) cancels those from {nH v + n"P) /{n-l) in <fy® . The other 
two terms g^ v and v?lH v j{n ■ I) 2 lead to gauge-independent results of order 
unity, which are thus neglected. 

Based on the above analysis, we need to compute only the contributions 
to from and from n 2 /^/(n • I) 2 . Details of the calculation refer to 
the Appendix. The anomalous dimension 

A«£ = -^C F (59) 

7T 

is then obtained, which is independent of n. The evolution of 4> is given by 



<P(^b,fi) = exp 



\^a s (yM)) 
i/(6M) y 



0(e,6,l/6) - (60) 



H{£,M,n) = exp 



Combining the evolution of the hard scattering to the largest scale M, 

f —\j > (a s (yM))-f 1 —\ q (a s (y/yM)) 
Jti/M y J^/m 2 y 

xH(C,M,M), (61) 

we derive 

^ = Jdtj d 2 h^ b, l/b)J(p', b, y/M/b)U(b, l/b)H(M, M) 

xexp[-S 7 (M,6)] , (62) 
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with the Sudakov exponent 



S 'y 



r 1 dy 

h/(bM) y 



pyM 2 rfn 2 

-^Xxiasifl)) + X^(a s (yM)) + X q (a s (^yM)) 

M 2 [i z 



' -(63) 



Note that the above expression differs from that employed in the PQCD 
analysis of inclusive B meson decays ||, where the lower and upper bounds 
of the variable fi are set to 1/b and yM, respectively. The distinction is 
attributed to the different ways to regulate the infrared divergences in K'\ 
The function 5((p' — l) 2 ) was introduced in the present work, while the Fourier 
factor e dT ' b was introduced in ||, such that the typical scales of K' are p'~ 
and 1/b, respectively. 

For comparision, we exhibit the result from the Wilson- loop approach [ffl, 



S 



N 



dy 

n /N y 



v m2 dkj 

y 2 M 2 k t 



2 T cusp {a s {k t )) + T{a s {yM)) + j(a s (^yM)) 

(64) 

where Sn describes the evolution of the iV-th moment of the differential 
decay rate, with the b quark mass M, the constant n = e~ 7B , 7^ being 
the Euler constant. The anomalous dimensions are given by r cusp = Xk, 
T = — (a s /n)CF and 7 = —a s /ir 0, ^3|. From Eqs. (|63| ) and (|64|), we easily 
identify the equality of the anomalous dimensions: 

= r cusp , X <j) = T , X q = 7 . (65) 

The two exponents S 7 and Sn are basically the same except for the lower 
bounds of y. The scale 1/b comes from the inclusion of transverse degrees 
of freedom in our formalism, and the N dependence is due to the moment 
analysis in H. 



VI. CONCLUSION 

In this paper we have discussed the resummation of large radiative cor- 
rections in inclusive processes at the end points, such as deeply inelastic 
scattering, Drell-Yan production and inclusive B meson decays, which lead 
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to the same results as those from the conventional Wilson-loop formalism. 
We have demonstrated the equivalence of the two very different approaches 
by carefully choosing the infrared cutoffs for the jet functions in deeply in- 
elastic scattering and in inclusive radiative B meson decays. Furthermore, 
we have shown in details that the resummation technique is gauge invari- 
ant: It can be performed both in covariant gauge and in axial gauge with 
an arbitrary gauge vector. The gauge dependent factors in the distribution 
functions and in the jet functions are canceled by those in the soft functions. 

We emphasize that a complete gauge-invariant set of higher-twist dia- 
grams contains those with valence partons carrying transverse momenta as 
considered in this work, and those with more partons entering a hard scatter- 
ing. We observe that the sole inclusion of the former in the study of inclusive 
processes at the end points is sufficient to guarantee the gauge invariance. 
Take the distribution function in deeply inelastic scattering as an example. 
The fractional parton momentum £Q, which is the upper bound of the vari- 
able p + in the resummation integral, approaches the largest scale Q at the 
end point. The upper bound of the scale [i in the evolution integral for the 
hard scattering is also Q, such that the cancellation of the gauge dependent 
factors is exact. 

In the region away from the end points there then exist residual gauge 
dependent factors in the form 



In order to maintain the gauge invariance, at least the gauge-dependent con- 
tributions from multiparton diagrams should be able to be resummed into 
an exponential, such that Eq. (|66|) can be canceled. This observation hints 
that it is possible to develop a new resummation technique for higher-twist 
effects from multiparton diagrams. We shall discuss this issue in a separate 
work. 

An immediate extension of the formalism presented here is to organize 
large corrections in processes involving gluons that enter a hard scattering 
as partons. This extension will be very useful for the study of, say, direct- 
photon processes, photoproduction and gluon fusion. We shall investigate 
this subject elsewhere. 




(66) 
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APPENDIX 



In this appendix we present some details of the one-loop evaluation of the 
functions K, G, K', G' and the anomalous dimension \$. 
1. The function K in Eq. ( |15| ) 

The function K from the two diagrams in Fig. 9(c) is written as 

K = ig C F fj. J ^ + 2vrz5(/ )e j ^ : - SK . 



v ■ nn ■ Iv ■ I 



Consider the former part of the integral without the exponential factor, 

d 4 ~H n 2 v^v u N^ 



K x = ig 2 C F [i e J 



(2vr) 4 - e vnn-lv IP 



The first term g^ v in N^ u does not contribute, since it leads to a numerator 
v 2 = 0. The second and third terms, —n^l u /(n ■ I) and —n u l^/(n ■ I), give 
equal contributions, which are given by 



K u = K 13 = -ig 2 C F ^n 2 J 



(2tt) 4 - £ P(n-iy 



We evaluate the above formula as a contour integral, assuming n + < 0, 
n~ > and riy = for simplicity. Certainly, the final expression should be 
Lorentz invariant. The reason for choosing a space-like n will become clear 
in the evaluation of G. Performing the integration over / + , whose singularity 
is pinched for l~ < 0, we obtain 

d 2 ~H r° „ 21- 
dl 



(2tt) 3 - £ J-oo {2n+l- 2 + n-l 2 ) 2 
Doing the l~ integration straightforwardly, we arrive at 

d 2 ~n 1 



K 



12 



-g 2 C F /j € J 



(2tt)3-^ l 2 T + 



in which the small constant a serves as an infrared cutoff, and will be set to 
zero at last. Computing the above integral using the standard dimensional 
regularization, we derive 

a s (2 V\ 
K12 = -— C F - + In - T — . 

2-n \e a 2 e^ E 
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The fourth term in , n 2 \H v jin ■ I) 2 , leads to the integral 

^ 2 ' d 4 ~n vi 



v ■ n J 

■ 2/7 

= 19 C F fx 



vn J (27r) 4 - e / 2 (n-/) 3 
n 2 ) 2 ( 1 d \ f d 4 ~ e l 1 



vn \ 2 dn- J J (2vr) 4 - e l 2 (n ■ I) 2 ' 

where the integral to be evaluated is the same as K i2 . It is then trivial to 
show = —K 12 , and thus K 1 = K 12 + K i3 + K u = K 12 . The second 
expression in Eq. @ also follows this conclusion. 

We now consider the latter part of the integral with the exponential factor, 

d 4 ~ e l -,. 2 ^i T . b "V^ 



K 2 = -g 2 C F ^J—^6(l 2 )e 



(2ir) 3 e v ■ nn ■ Iv ■ I 

Following the above procedures, we find that only the contribution from the 
second term of needs to be computed, which gives a similar expression 
to K 12 , 

r d 2 ~ e l 1 
K2 = 9 C ^ J ( 27 r)3-^ + fl 2 e • 

Since this integral is free of ultraviolet and infrared singularities, it can be 
evaluated by setting e = directly. We obtain 

K 2 = ^C F K (ab) « -^C F \n—^ 

in the limit a — > 0, where K is the modified Bessel function of zeroth order. 
The sum K\ + K 2 is then given by 



K l + K 2 = -^C F 



- + ln(7r / u 2 6 2 e 7E ) 
e 



which is independent of the cutoff a. If renormalizing K in MS scheme, ie., 
choosing 



5K = _^ Cf ( 2 - 

2tt Ve 



+ In 47T - 7 B j , 
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we have the lowest-order K = K\ + Ki — 8K as 



K = -^C F 

71 



ln(6/i) + In ■ 



Neglecting the second term of order unity, we derive Eq. (JT 



2. The function G in Eq. ([16]) 

The function G from Fig. 9(d) is written as 



G = ig 2 CpH e J 



d 



4-e 



I ■nfVpNi" 



(27r) 4 " e v-nn- IP 



(p+iy 



V-l 



5G . 



A simple investigation shows that only the third term —n u l^/ (n ■ I) con- 
tributes. Hence, we express the integral as G13 — K13, where Ki 3 = K u has 
been computed before, and G13 is the former part of the integral, 

i2 " d A ~n {t>+ jt)vi 



G 



13 



-ig C F /j e 



n 



v ■ n 

3-n ..e„2 



-2ig l C F n e n l 



(27r) 4 - e (p + l) 2 (n ■ Ifl 2 ^ 
d 4 ~H (p + i)+r 
(2tt) 4 - £ (p + l) 2 (n-l) 2 l 2 ' 



In the second expression those terms, vanishing after contracted with the 
matrix structure of the initial hadron, have been neglected. Performing the 
integration over l~ and then over Z + , we obtain 



G 



13 



g 2 C F tf 



n 



d 2 ~n 



2p+ J (27r) 3 - e 
2n + p ■ n 



n + l 2 T 



■In 



n 2 [(p ■ n) 



n 2 ll' 



(p ■ n) 2 
2p+ 
n 2 {l 2 T + a 2 ) 



n 2 l 2 



J (p ■ n) 2 — n 2 l\ + p ■ n 



n 2 l 2 



p ■ n 



It is obvious that choosing n 2 < avoids extra singularities in the above 
integral. The first and second terms have no infrared divergences, but the 
third term does. Hence, we have introduced a small cutoff a into the de- 
nominator of the third term as in the evaluation of Ki 2 . After a tedious 
calculation, we arrive at 



G 



13 



2tt a 2 e 



2 v 2 
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which is ultravioletly finite because of the exact cancellation of poles. At 
last, G13 — K13 is written as 



G 



13 



K 



13 



— + ln^ - 

e 7r/i z 



Again, the dependence on the cutoff a disappears. In the MS scheme, we 
have 5G = —SK and G = G13 — K13 — 5G as 



G = -^C 

71 



In h In 



Keeping only the first logarithmic term, we derive Eq. (|T8|). 



3. The function K' in Eq. 

Based on the conclusion for K that only the contribution from the term 
nH v I (n ■ I) needs to be evaluated, we reexpress Eq. (|29| ) as 



5(p' 2 )K' = -ig 2 C F ^n 2 J 
-8(p' 2 )8K' . 



d A ~ e l 



(2vr) 4 - e (n • If 



5{p 



12 \ 



2ni5(l 2 )5((p' - If 



Since the first term has been computed, we concentrate on the second term. 
Similarly, the integration over l + leads to 



5{p' 2 )K' 22 = g 2 C F ^n 2 J 



d 2 ~H 
(27r) 3 - £ 



21 



n-l 2 T ) 2 



The integration over I gives 
5(p' 2 )K' 22 = g 2 C F ^n 2 J 



d 2 ~n 



{2n+l- 2 -f 

2{p'-) 2 p' 2 



8W- 



1- 

p F. 



(2vr) 3 - e [2n + (p'-) 2 (l 2 T + a 2 ) + n-p' 4 } 2 ' 



where the small constant a has been inserted for the similar reason. Because 
the above integral does not contain ultraviolet and infrared singularities, we 
set e = and perform the integration over It- We obtain 



a. 



5(p' 2 )K' 22 = -C F 



V 



12 



n p'* + 4is' 2 (p'-) 2 a 2 
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Now we interpret the right-hand side of the above formula as a distribu- 
tion function, and extract the on-shell pole at p' 2 = using 

m r 2 (p'~) 2 n' 2 



5(p' 2 )K> 2 = ^C F 5(p' 2 ) f " dp 

It JO 



p' 4 + 4z/ ,2 (p'") 2 a 2 



= 6(p' 2 )^C F In ^ . 
Combining the above expression with Ku, we derive Eq. (^1 

ix p' 



4. The function G' in Eq. (H) 

The calculation of G' from Fig. 10(d) is basically the same as of G. The 
treatment of the third and fourth diagrams has been explained in Sec. III. 
That is, we replace p'~ by p + and a by y / p + p'^ in —K' 22 . 



5. The anomalous dimension A^ in Eq. fl59|) 

We evaluate the anomalous dimension A^ of the distribution function 
in inclusive B meson decays. As stated in Sec. V, single soft logarithms are 
produced from Figs. 6(b) and 6(c). To extract the soft structure of these 
radiative corrections, we apply the eikonal approximation to the b quark line 



lOfl . Fig. 6(b), corresponding to the self-energy correction to the b quark, 



gives a contribution proportional to 



d 4 l 



(2tt) 4 v ■ ll 2 



where 7^ comes from the outer gluon vertex on the b quark line. This integral 
vanishes because the integrand has an odd power of I. 

The contribution from Fig. 6(c) has been expressed in Eq. (0). As ex- 
plained in Sec. V, only the terms g^ v and n 2 l^l u / (n-l) 2 need to be considered, 
which lead to equal results. Hence, we have the bare distribution function 

d 4 ~H n 2 
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= 2g 2 C F ^J 
from which the infrared pole 



d 2 ~ e l 1 .. 



(2vr) 3 - e l\ 



e~ Th 



AO) _ a v 1 

is extracted. Therefore, the anomalous dimesion is given by 



7T 

as exhibited in Eq. (|59|). 



= Cf , 
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Figure Captions 

Fig. 1. (a) The jet subprocess defined in Eq. P). (b) and (c) Lowest-order 
diagrams of (a). 

Fig. 2. The derivative p + dJ/dp + in covariant gauge. 

Fig. 3. (a) 0{a 2 s ) examples for the differentiated J. (b) Factorization of 
the function K at 0(a s ). (c) Factorization of the function K at 0(a 2 s ). (d) 
Factorization of the function G at 0(a s ). 

Fig. 4. (a) General diagrams for the functions K and G. (b) The 0(a s ) 
function K. (c) The 0(a s ) function G. 

Fig. 5. (a) The derivative p + dJ/dp + in axial gauge, (b) General diagrams 
for the functions K and G. (c) The 0(a s ) function K. (d) The 0(a s ) 
function G. 

Fig. 6. (a) Lowest-order diagram of deeply inelastic scattering, (b), (c), (d) 
and (e) 0(a s ) corrections to (a). 

Fig. 7. Factorization of Fig. 6(e) for (a) the loop momentum parallel to the 
momentum of the valence quark, for (b) the loop momentum parallel to the 
momentum of the outgoing quark, and for (c) the soft loop momentum. 

Fig. 8. Factorization of deeply inelastic scattering. 

Fig. 9. (a) The derivative p + d<p/dp + in axial gauge, (b) General diagrams 
for the functions K and G. (c) The 0(a s ) function K. (d) The 0(a s ) 
function G. 

Fig. 10. (a) The derivative p'~ dJ / dp'~ in axial gauge, (b) General diagrams 
for the functions K' and G '. (c) The 0(a s ) function K' . (d) The 0(a s ) 
function G' . 

Fig. 11. (a) Lowest-order diagram of Drell-Yan production, (b), (c), (d) 
and (e) 0(a s ) corrections to (a). 

Fig. 12. Factorization of Drell-Yan production. 
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